The Hyperfine Splittings in Bottomonium and the B q (q = n,s,c) Mesons 
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A universal description of the hyperfine splittings (HFS) in bottomonium and the B q (q — n, s, c) 
mesons is obtained with a universal strong coupling constant a s (fi) = 0.305(2) in a spin-spin po- 
tential. Other characteristics are calculated within the Field Correlator Method, taking the freez- 
ing value of the strong coupling independent of n/. The HFS M(B*) - M(B) = 45.3(3) MeV, 
M(B*) — M(B 3 ) — 46.5(3) MeV are obtained in full agreement with experiment both for nj = 3 
and n/ = 4. In bottomonium, M(T(9460)) - M(r/ b ) = 70.0(4) MeV for n f = 5 agrees with 
the BaBar data, while a smaller HFS, equal to 64(1) MeV, is obtained for n/ = 4. We predict 
HFS M(T(2S)) - M{ Vb (2S)) = 36(1) MeV, M(T(3S)) - M(r/(3S)) = 27(1) MeV, and M(B*) - 
M(B C ) = 57.5(10) MeV, which gives M(B*) = 6334(1) MeV, M(B C (2 1 S )) = 6865(5) MeV, and 
M(B*{2S 3 S 1 )) = 6901(5) MeV. 
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I. INTRODUCTION 

Recently, r]b(lS) has been discovered by the BaBar col- 
laboration in the radiative decays T(3S) —> jrj b (lS) |l| 
and T(25) — > ■7775(15') [2j, with a mass (averaged over 
two results) M (r) b ) = 9391.1 ± 3.1 MeV. It gives a rather 
large hyperfine splitting (HFS), A(bb) = M(T(1S)) - 
M(rjb(lS)) = 69.9 ± 3.1 MeV. Later this mass was con- 
firmed by the CLEO collaboration also in the radiative 
T(3S) -> jr) b (IS) decay |[. This important new in- 
formation allows to test again our understanding of the 
hyperfine (HF) interaction in QCD. 

Although a spin-spin potential between heavy quarks 
was used in numerous studies, the parameters defining 
this potential significantly differ in different models. As 
a result, theoretical predictions for the mass difference 
A(bb) = M(T(9460)) - M(r) b (lS)) vary in a wide range: 
35 — 90 MeV Q-Q and in most cases they are smaller 
than the experimental number. On a fundamental level 
the spin-spin potential V ss (lat) has been studied in detail 
in quenched QCD (see and references therein). This 
lattice potential appears to be compatible with zero at 
distances r > 0.30 fm and (for unknown reasons) has neg- 
ative sign at smaller r (with a large magnitude); in any 
case the lattice potential does not contradict the Fermi- 
Breit potential with S 3 (r) although the behavior of 
the spin-spin potential at r < 0.3 fm remains uncertain. 

On the other hand, a detailed phenomenological anal- 
ysis given in Ref. 4] has demonstrated the importance of 
the smearing of the <5 3 (r)-function, from which one may 
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expect that for heavy mesons, containing a 6-quark, the 
use of 5 3 (r) may be a good approximation. For lighter 
mesons, like D, D s , and charmonium, a nonperturbative 
spin-spin potential may be essential, giving a contribu- 
tion ~ 10% [HI]. Here we concentrate on bottomonium 
and the B q (q = n,s,c) mesons, for which nonpertur- 
bative contributions are small, and neglect the smearing 
effect, in this way avoiding to introduce several unknown 
parameters. 

Our main goal here is the extraction of the strong cou- 
pling cvhf (n) from known HFS. In theoretical models two 
typical choices of evhf are used: 

1. In the first one, "a universal" ahf is used. For exam- 
ple, in Ref. [j| a^f = 0.36 taken, was obtained from 
a fit to the mass difference M(J/ip) - M (f/ c (15)) 

1 1 7 MeV, but their HFS, Af(T(9460)) - M(r? 6 ) = 
87 MeV, is ~ 25% larger than the experimental 
number. In Ref. @, using a smaller a^s = 0.339 
a good description of the HFS of the B and B s 
mesons was obtained. However, a comparison of 
their and our results is difficult, because a large 
string tension, a = 0.257 GeV 2 , was taken in Q, 
while here and in Ref. [4| the conventional value 
a = 0.18 GeV 2 is used. 

2. The second choice, with a scale \i dependent on 
the quark mass, is mostly used in pQCD, where 
Q!hf(wh) ~ 0.18 and evhf(w c ) ~ 0.26. Just due to 
such a small value of a n f(™&), taken in bottomo- 
nium, small HFS were obtained in Ref. [T2J], al- 
though their wave functions (w.f.) at the origin 
gave excellent descriptions of the dielectron widths 
for T(nS)(n= 1,2,3) El- 
Here we use instead of the Fermi-Breit potential a spin- 
spin potential derived using the Field Correlator Method 



2 



(FCM) [lj], where relativistic corrections are taken into 
account and with the mass of a light quark m n = 5 MeV 
(n = u, d) and m s = 200 MeV for an s-quark, the B, B s 
mesons can be considered on the same footing as the B c 
mesons. 

It can be shown that the HFS are sensitive to the value 
of A^jg(n/) taken. Since A^jg is known only for rif — 5 
and Ajyfg, used for 77/ = 3, 4, varies in a wide range, we 
make here the assumption, already used in that the 
freezing value of the vector coupling constant (denoted 
as a cr it(rif)) is the same for 77/ = 3, 4, 5. 

Then it appears possible to obtain a good description 
of the HFS for the B q mesons (q — n,s,c) and bottomo- 
nium, taking a universal a s (/i) = 0.305(2) (with small 
one- loop corections); its value is smaller than in Refs. [5j 
and [9|. 

We also predict the HFS of the as yet undiscovered 
?7t, (25") and r/b(3S), and the mass of B*. 



II. THE HF POTENTIAL IN THE FIELD 
CORRELATOR METHOD 



The Fermi-Breit potential, 

32tt a hf (ju) 



V ss (r) = si • s 2 —— — 



777,17712 



(1) 



widely used in heavy quarkonia, contains the constituent 
quark masses rh\ and rh-x-, which are very much model- 
dependent. 

In Eq. (flj the strong coupling constant a^fOu) may 
differ from a s (/_t) (in the MS renormalization scheme) 
due to higher order perturbative corrections. Here we 
take the one-loop corrections from [TBI ]. 



Ohf(#) = a s (fi) 
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where 



5 n 8 3, „ 

p = — B ln2. 

y 12 H0 3 4 



(2) 



(3) 



Although higher corrections are small: ~ 6% for 77/ = 
3, - 3% for n/ = 4, and < 0.1% for 77/ = 5, still they 
are needed to improve the accuracy and have a more self- 
consistent picture. 

The important role of relativistic corrections, even for 
the B c meson, has already been underlined in Refs. [4] 
and Mi, a nd recently in the lattice calculations of the B* 
mass [lj| , [I?} • We take them into account using the spin- 
spin potential (without smearing), derived in the FCM 
G3, El: 



T> I \ 32?r a hf(M) c3/ \ 

V ss {r) = si • s 2 — <r(r), 
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for which the HFS is 



A hi (nS) = 8 -^\R n (0)\ 2 . 

9 UJ\U)2 



(4) 



(5) 



In Eqs. (HJ and ([5]) the variables uji(nS),u)2(nS) are the 
averaged kinetic energies of a quark 1 and a antiquark 2, 
which play a role of the dynamical masses: 



(Ji(nS) = (\f p 2 + ml) nS , W2(nS) = {\J p 2 + rn|)„s. 

The important point is that in Eq. (|6|) the masses mi and 
77i2 arc well defined; they arc the pole masses of c and b 
quarks (known now with an accuracy of ~ 70 MeV for a 
b quark and ~ 100 MeV for a c quark [lj|). In leading 
order, the pole mass does not depend on the number of 
flavors, while to order (a s (?77Q) 2 ) it slightly depends on 
77/ . We take here mi = m n = 5 MeV for a light quark 
(71 = u,d); 777 s = 200 MeV for an s quark; the pole 
mass TO C = 1.41 GeV; nib = 4.79 GeV for 77/ = 3 and 
TO,, = 4.82 GeV for 77/ = 4, 5. 

The quantities uji and the w.f. are calculated with 
the use of the relativistic string Hamiltonian (RSH) , also 
derived in the FCM [13], 



P 



where the variables Ui enter as the kinetic energy op- 
erators. However, if one uses an einbein approximation 
[l~8j . f2l"j and considers the spin-dependent potential as 
a perturbation, then oji should be replaced by its matrix 
elements (m.e.) ((6|). 

A simple expression for the spin-averaged mass M(nS) 
follows from the RSH [U: 

777? uio rnk „ 
M(ns) = ^ + ^- + ^ + ^-+E nsM . (8) 

Here, the excitation energy E n s(uj rc ^) depends on the 
reduced mass: uj- 
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The mass formula (151) 
does not contain any additive constant in the case of 
bottomonium, while for the B and B s mesons a nega- 
tive (not small) self-energy term, proportional to (ui q )~ l 
(q = 77, s), has to be added to their masses [22l |. 

Then the variables uii(nS), the excitation energy 
Ensured), and the w.f. are calculated from the Hamilto- 
nian (JZ} and two extremum conditions, d M(nS) 1 1 dujj = 
(i = 1, 2), which are put on the mass M(nS) fl8l |: 



2 2 2wi 

= E(nS)(f nS , 

L0 2 (jlS) = 777 2 



2cj 2 



2w r 
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V B (r) 



<Pns(r) 



dE(nS, /ire, 
du>i(nS) 



-, (7 = 1,2). 



(9) 



(10) 



In a Hamiltonian approach the choice of the static po- 
tential Vb(t) is of great importance; we take it as a sum 
of a linear confining term and the OGE -type term: this 
property of additivity is well established now in analyti- 
cal studies and on the lattice [23l|. [24]: 



V B {r) 



a r ■ 



4obW 
3 r ' 



(11) 
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For the string tension a conventional value, a = 
0.18 GeV 2 , is used here for all mesons. 

The main uncertainty comes from the vector coupling 
ay(r), which is taken here from Refs. [25| . (26j j and de- 
noted as aB(r). Two important conditions have to be 
put on the vector coupling: 

(i) As in pQCD, it must possess the property of asymp- 
totic freedom (AF); precisely owing to this property the 
static interaction depends on the number of flavors. 

(ii) The vector coupling freezes at large distances. The 
jropert y o f freezing was widely used in phenomenology 
1-0, H3 and ol 
static potential [25 



propert y q t 

|3|-0) |23| and observed in lattice calculations of the 



Unfortunately, one cannot use the static potential and 
the freezing (critical) constant from lattice studies, where 
the latter is found to be significantly smaller than in phe- 
nomenology and background perturbation theory (BPT) . 
There a B (crit) = 0.58 — 0.60 is used (these numbers are 
close to the value from Q). On the lattice, a cr it(lat) ~ 
0.30 in full QCD (rif = 3) J2I and a crit (lat) - 0.22 in 
quenched calculations (loj . |23| were obtained. 

In Eq. (|lip the vector coupling «b (t) is defined via the 
vector coupling a B (g 2 ) in momentum space (26j : 



2 f sm(gr) 

it J q 





(12) 



where the vector coupling a B (q 2 ) is taken in two-loop 
approximation, 



4tt / fix lnt B 



PotB \ Pa *B 



with the logarithm 



Ml 



A 2 



(13) 



(14) 



containing the vector constant Ab (rif), which differs from 
the QCD constant A^m(rif). The relation between them 
has been established in Ref. 1281: 



A B (n/) = Aj^exp 



ai 
Wo 

10 , 



(15) 



with /3q = 11 — |n/ and a\ = # — -^n/. Therefore the 
constant Ab is always larger than A^g, 



A 



(5) 



1.3656 Af^n, 



5), 



V„" - 1.4238 Ai^nj =4), 



^(4) 
V B 

,(3) 



(16) 



A',;" - 1.4753 A^n/ = 3). 



At present, only the QCD constant Aj^g(n/ = 5) is 
known with a good accuracy, while for rif = 3,4 it 
is defined with an accuracy ~ 10% (l9j). For a given 



TABLE I: The kinetic energies ui q (lS) (q = n, s, c) and oii,(15) 
(in MeV) for the static potential Vb(?") with n/ = 4 and 
Qcrit =0.58 



Meson S 


B a 


B c 


TOq 5 


200 


1410 


CJq — r/lq 611 


489 


248 


Wf, 4805 


4825 


4888 


Ub — mb 25 


25 


83 



TABLE II: The ratios g Bq JTUJ (in GeV) and (0) | 2 



(in 



for the B,(1S) 


mesons 


( Qcrit = 


0.58,71/ = 4) 




B 


B s 


B c 


|i?(0)| 2 


0.477 


0.551 


1.669 


5s„ 


0.161 


0.166 


0.205 



^Ms( n f ~ 5) one can define a C rit("-/ = 5). Then, assum- 
ing that the freezing constant is the same for rif — 3,4, 
the QCD constant, A^jg for rif — 3, 4, can be calculated. 

The mass Mb under the logarithm £b is proportional 
to yfa and for a = 0.18 GeV 2 its value is equal to Mb = 
1.0 ±0.05 GeV H|. 

The critical constants «B(crit) = a B (q 2 = 0) in mo- 
mentum space and an(r — > 00) in coordinate space are 
proved to be equal [261 ] . 

Our calculations give small relativistic corrections for 
bottomonium: u> b (lS) - m b ~ 190 MeV (~ 4%) and 
~ 7% for the 25* and 3S states. It is of interest to notice 
that the relativistic correction to the 6— quark mass is 
even smaller in the B c meson, uj b (lS) — m b ~ 83 MeV 
(~ 2%), while for a c— quark the difference oj c (1S) — m c ~ 
250 MeV is already - 20%. The values of u q (lS) - m q 
are given in Table U for the B q mesons (q — n,s,c). 

It is convenient to introduce the ratio gs , 



9B q (nS) 



K(o)l 2 



ui\ (nS)ui2(nS) ' 



(17) 



which directly enters the HFS (O and appears to be 
weakly dependent on small variations of the masses mi 
and mi , which are compatible with a good description of 
the meson spectrum. 

The w.f. at the origin are sensitive to the values of 
Ab(»v) ( m only the case with rif = 3 was consid- 
ered). However, if the same freezing value of the coupling 
constant is taken for rif — 3, 4, 5, then the differences be- 
tween the w.f. at the origin for rif — 3 and rif = 4 are 
< 3%. 

We use here two values for a cr ;t: a C rit = 0.580 
and 0.604, for which corresponding values of Ab(«/), 
Ajyjs^/) are given below in Eq. (TiUl) . 

In bottomonium the difference between g b for rif = 4 
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TABLE III: The ratios g b {nS) (in GeV) and |i?„(0)| 2 (in 
GeV 3 ) for the 15, 2S, and 3S bottomonium states with 
acrit = 0.604 and for n/ = 4, 5 





IS 


2S 


3S 


l^(0)| 2 (n/ = 


5) 6.476 


3.398 () 


2.682() 


g b {n f = 5) 


0.258 


0.134 (1) 


0.105(1) 


\R n (0)\ 2 (n f = 


4) 5.668 


3.126 


2.508 


g b (n f = 4) 


0.230 


0.127 


0.100 


A, l/ (n / = 5) 


70.3 


36.3 


28.4 



and rif = 5 appears to be larger, ~ 10% (see Table HTH) . 
where in both cases a cr it = 0.604 is used. 

For the values of a cr it(n/) = 0.604(0.58) and 
A m (n f = 5) = 0.245(0.236) GeV, the following vec- 
tor constants: A B (n/ = 3) = 0.40(0.389) GeV; A B (n f = 
4) = 0.372(0.360) GeV, A B (n f = 5) = 0.335(0.323) GeV 
are obtained. Then from the relation (fT5|) we have 

a ms("/ = 3) = 271(264) McV, 

A M g(n / = 4) = 261(253) MeV, (18) 

X m {n f = 5) = 245(236) MeV. 

For rif = 5 it gives a s (Mz) (two — loop) = 0.1194 for 
a crit = 0.604 and a s (M z ) = 0.1188 for a crit = 0.58; both 
numbers agree with the world averaged value, a s (Mz) = 
0.1176 ± 0.0020 within its error bar. 

III. RESULTS 



A(66) is equal to 36(2) MeV and 27(2) MeV, respectively, 
for the 3S and 2S states. 

For the B q mesons, the cases with rif = 3, 4 have been 
considered and in both cases agreement with experiment 
is obtained (see Table llV)) : it happens because the by 
~ 3% smaller value of gB q is compensated by the ~ 3% 
larger value of (Xhf (n>f = 3). Thus for B and B s the 
number of flavors cannot be fixed, if only data on the HFS 
(and the spectrum) are fitted; therefore some additional 
information, like decay constants, is neededd to fix n/. 



TABLE IV: The HFS (in MeV) of the B q mesons with 
a s (fi) = 0.307 (a h f(n/ = 4) = 0.316, a M (n f = 3) = 0.324) 





B 


B s 


B C (1S) B C (2S) 


A h f(n/ = 


4) 45.2 


46.5 


57.8 38.0 


Ahf(n/ = 


3) 45.4 


46.1 


57.2 37.4 


Ahf(exp) 


45.78 ± 0.35 


46.5 ± 1.25 


absent absent 



From our calculations we predict the following masses 
of the triplet and singlet c6(2S) states: M(B*(2 3 Sx)) = 
6901(5) MeV and M(B C {2 x 5 )) = 6865(5) MeV, which 
are calculated in single-channel approximation. 

Notice that the extracted strong coupling constant, 
ct s (/j,) — 0.305(2), is smaller than the one used in Refs. @ 
and [9j. The renormalization scale, /x ~ 1.6 — 1.65 GeV, 
corresponding to this coupling, is rather large but agrees 
with the existing interpretation of the spin-spin potential 
as dominantly perturbative, thus partly justifying the use 
of the <5 3 (r)-function. 



The experimental error in the HFS 

A hf (66) = M(T(946Q)) -M{r, b ) = 69.9 ±3.1 MeV (19) 

is small, ±3 MeV, and it is even smaller, < 1 MeV, for 
the mass differences M(B*) - M(B), M(B* S ) - M(B B ) 
[T^ |. Therefore, we expect that the coupling constant ckm 
can be extracted with a good accuracy from these data. 

Notice that from Eq. |(3J) and the value a s (fi) = 
0.305(2) one obtains the following coupling constants 
with one- loop corrections: a^i(nf = 5) = a s (ii) — 
0.305(2), a hf (n f = 4) = 0.314(2), and a hl (n f = 3) = 
0.323(2). Precisely these values are used in our analysis. 

For g b = 0.230 GeV (n f = 4) and g b = 0.258 GeV 
(rif — 5) (see Table llllj) and taking Qm(^/ = 5) = 
a s (fi) = 0.305(2) and a M (n f = 4) = 0.314(2), we ob- 
tain Am(66) = 64.2(4) MeV (n f =4) and 70.0(4) MeV 
for nf — 5. The difference between them, ~ 10%, is not 
small and one may conclude that the HFS in bottomo- 
nium is in full agreement with the BaBar data [l|, @ 
only for n/ = 5. 

For the 25* and 3S bottomonium states the difference 
between the cases with rif = 4 and rif = 5 is small; their 
HFS coincide within 2 MeV. Therefore we predict that 



IV. CONCLUSION 

Our study of the HFS is performed assuming that the 
freezing value of the coupling constant is the same for 
rif = 3,4,5 and considering a cr i t = 0.58 and 0.60. The 
calculated HFS of the B and B s mesons are in good 
agreement with experiment for both freezing constants. 
It happens that the HFS for rif — 3 and nj — 4 co- 
incide with each other (within 0.5 MeV), if the one- 
loop correction is taken into account in ahf(n/). The 
HFS, averaged over two results, are M(B*) — M(B) = 
45.3(2) MeV, M(B* S ) - M(B S ) = 46.5(3) MeV, and 
M(B*) - M(B C ) = 57.5(10) MeV. The latter num- 
ber gives the mass of the as yet unobserved B* me- 
son, M(B*) = 6.334 ± 1 MeV. For excited B C (2S) states 
we predict the masses: M(B*(2S) = 6901(5) MeV and 
M(B C (2S) = 6865(5) MeV, which are calculated neglect- 
ing open channels. 

In bottomonium the choice rif = 5, where Am (66) = 
70.0(4) MeV, is considered preferable, since for nj = 4 
the HFS is smaller, equal to 64.2(4) MeV . 

For the 2S and 3S bottomonium states, the calculated 
HFS are 36(1) MeV and 27(1) MeV, respectively. 
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